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ABSTRACT 


The  final  phase  in  the  development  of  a  low  altitude  turbulence 
model  for  ai*.  craft  gust  load  application  is  described.  The  basic  data  used 
for  this  model  are  the  B-66B  low-level  gust  data  and  climatological  wind 
data  provided  by  the  National  Weather  Records  Center,  Asheville,  N.  C. 

In  preceding  studies  the  turbulence  spectrum  functions,  rms  gust  velocity/ 
mean  wind  speed  functions,  and  mean  wind  speed  distribution  functions 
were  determined  for  variations  in  height,  atmospheric  stability,  and  ter¬ 
rain  roughness  conditions,  hi  this  final  study  needed  to  complete  the  tur¬ 
bulence  model,  climatological  wind  statistics  are  used  to  relate  average 
mean  wind  speed  characteristics  to  terrain,  height,  thermal  stability,  time 
of  day  and  seasonal  variations.  The  data  represent  locations  throughout 
the  United  States.  Based  on  results  obtained  from  the  climatological  wind 
statistics,  procedures  are  outlined  for  applying  the  turbulence  model  to 
estimate  aircraft  gust  load  exceedances  for  a  specified  low  altitude  opera¬ 
tional  history.  A  preliminary  comparison  of  the  model  and  B-5Z  service 
load  data  is  made  for  an  800  ft  terrain-clearance  altitude. 


This  abstract  is  subject  to  special  export  controls  and  each  trans¬ 
mittal  to  foreign  governments  or  foreign  nationals  may  be  made  only  with 
prior  approval  of  Air  Force  Flight  Dynamics  Laboratory,  FDTR,  Wright- 
Patterson  Air  Force  Base,  Ohio  45433. 
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NOMENCLATURE 


tp  normalized  spectrum  function 
f  frequency 

k  wave  number,  f/U 

O  wave  number  times  Zv 

T(ifl)  transfer  function  of  any  airplane  response  item  to  sinusoidal 
gust  of  unit  amplitude 

U  mean  wind  speed  or  airplane  forward  velocity 

u  mean  wind  speed 

u  climatic-wind-speed 

h,  z  height  above  surface 

aw  rms  of  vertical  gust  velocity 

L  scale  parameter 

Y  general  load  parameter 

N(Y)  number  of  exceedances  of  a  given  load 

Nq  number  of  crossings  of  Y  of  zero  reference  value  with  a  positive 

slope  per  foot  distance  flown 

2 

A  integral  of  product  of  the  normalized  spectrum  function  and  the 
square  of  the  transfer  function 
H  hour  of  day 
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Section  1 


INTRODUCTION 

This  report  summarizes  the  final  phase  of  a  study  to  develop  a  tur¬ 
bulence  model  to  relate  low-level  atmospheric  turbulence  spectra  to  mean 
meteorological  parameters,  terrain  conditions,  and  height.  Previous 
phases  associated  with  the  model  development  are  reported  in  refs  1-4. 

The  model  is  designed  to  provide  a  gust  load  procedure  to  account  for  the 
effects  of  turbulence  on  aircraft  design  load  factors  and  for  estimating 
structure  fatigue  life.  The  procedure  requires  assigning  estimates  of 
diurnal  and  seasonal  flight  operations  (weights,  speeds,  etc)  in  conjunc¬ 
tion  with  terrain  and  altitude  conditions. 

A  summary  of  the  previous  three  phases  repor  ed  in  refs  1-3  fa  pre¬ 
sented  below: 

Ref  1:  Data  obtained  from  the  B-66  low-level  gust  study  (ref  5) 
were  analyzed  from  a  meteorological  and  a  power  spectral  density  stand¬ 
point.  Meteorologically,  the  basic  B-66  wind  speed  and  temperature  data 
were  combined  with  synoptic  observations  to  provide  wind  speed  and  lapse 
rate  estimates  consistent  with  both  B-66  estimates  and  synoptic  data.  From 
the  spectrum  data  analysis  a  spectral  shape  representation  with  height 
was  obtained  for  two  (rough  and  smooth)  terrain  classifications. 

Ref  2:  Relations  between  the  standard  deviation  of  vertical  velocity 
and  mean  wind  speed  for  smooth  terrain  conditions  were  obtained  for  three 
atmospheric  stability  classes;  and  1400ft  tower  data  were  analyzed  to  pro¬ 
vide  mean  wind  profile  characterisics  ' >r  the  three  lapse  rate  conditions. 
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Climatological  wind  statistics  for  Dayton,  Ohio  were  analysed  and  the 
possibility  of  obtaining  a  generalised  mean  wind  distribution  function  based 
on  the  average  climatic  wind  speed  was  suggested. 

Ref  3:  Climatological  data  for  smooth  and  rough  terrain  stations 
were  analysed  and  generalised  wind  distribution  functions  were  obtained. 
Procedures  for  estimating  the  rms  gust  velocity  probability  density  dis¬ 
tribution  functions  and  the  exceedance  probability  distributions  were  deter¬ 
mined. 

In  the  current  final  phase  of  the  basic  model  development,  exten¬ 
sive  climatological  data  are  analysed  for  the  purpose  of  (1)  confirming  the 
wind  speed  probability  density  distribution  functions*  and  (2)  providing  a 
climatic  mean  wind  speed  model  to  account  for  diurnal  and  seasonal  changes 
as  related  to  height,  atmospheric  stability,  and  terrain  conditions.  In  ad¬ 
dition,  step-by-step  procedures  are  outlined  for  applying  file  turbulence 
model  to  estimate  aircraft  gust  load  exceedance  probabilities. 
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Section  2 

GENERAL  APPROACH  TO  THE  PROBLEM 

A  summary  of  the  meteorological  and  geophysical  aspects  of  the 
low  altitude  gust  load  model  is  described  in  this  section. 

1.  Basic  load  expression 

The  basic  equation,  presented  in  refs  3  and  4,  describing  the 
average  number  of  peak  values  per  unit  distance  exceeding  a  given  load 
T  is  given  by 

*  2 

N(Y)  =  Nc  $  «dw)  exp  {  da„  (1) 


where  is  the  rms  of  vertical  gust  velocity,  f(ffw)  is  the  probability 

density  function  of  o  ,  and  N  and  are  defined  as 

wo 


To  determine  the  load  exceedance  probabilities  for  a  range  of 
aircraft  operating  conditions  {weight,  speed,  etc)  and  geophysical  en¬ 
vironmental  conditions  (wind  speeds,  atmospheric  stability,  topography, 
etc),  it  is  necessary  to  evaluate  Equation  (1)  for  significant  changes  in 
these  conditions.  Thus, 

N(Y)  =  Y  Ty  N(Y)  (2) 

where  ^  Ty  =  1  and  N(Y)  is  the  load  exceedance  probability  averaged 
over  conditions  ij. 

2.  Basic  approach 

The  basic  approach  to  the  low  level  gust  procedure  involves  model¬ 
ing  the  vertical  gust  velocity  intensity  (ff^)  and  spectrum  characteristics 

(up  )  in  terms  of  the  significant  geophysical  and  meteorological  condi- 
WN 

tions.  hi  general,  the  meteorological  conditions  include  mean  wind  speed 
and  atmospheric  stability  in  relation  to  topography,  height,  time  of  day, 
and  seasonal  parameters,  hi  the  previous  studies  (refs  1-4),  it  was  con¬ 
cluded  that  the  interdependence  of  some  of  these  conditions  was  sufficiently 
weak  to  permit  a  relatively  straightforward  application  of  Equation  (1). 

A  summary  of  these  conditions  is  as  follows: 

1)  The  form  of  the  normalized  vertical  gust  spectrum  (<p  )  is  considered 

WN 

independent  of  atmospheric  stability  when  the  frequencies  are  divided  by 
the  true  airspeed;  the  spectrum  form  depends  only  on  the  height  above 
the  surface  and  the  terrain  roughness  (see  refs  1  and  2). 


2)  Over  smooth  terrains,  the  rms  turbulence  intensity  (aw)  is  considered 
to  be  principally  related  to  mean  wind  speed  and  parametrically  related 
to  atmospheric  stability.  Over  rough  terrains  the  rms  gust  velocities 
are  considered  independent  of  stability  (refs  3  and  4). 

3)  The  probability  distributions  of  the  climatological  wind  speeds  are  con¬ 
sidered  independent  of  height  and  atmospheric  stability,  and  weakly  de¬ 
pendent  an  terrain  conditions,  when  the  mean  wind  distributions  are 
normalized  by  the  appropriate  average  mean  wind  speeds  representing 
height  and  stability  classifications  (see  refs  3  and  4).  This  wind  speed 
behavior  enables  the  f(ffw)  function  to  be  specified  on  the  basis  of  the 
rms  turbulence  intensity  function  and  the  average  mean  wind  speed. 

hi  Section  6  of  this  report,  the  average  mean  wind  or  climatic  wind 
speed  (as  it  is  referred  to  in  refs  2-4)  characteristics  are  considered  for 
representative  locations  in  the  United  States.  These  wind  data  are  sum¬ 
marized  for  both  smooth  and  rough  terrain  conditions  for  such  factors  as 
atmospheric  stability  and  height  above  the  surface  for  diurnal,  seasonal, 
and  yearly  periods. 
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Section  3 


SPECTRUM  SHAPE 

The  form  of  the  normalised  spectrum  function  is  described  in  refs 
3  and  4  for  the  low-level  regime.  The  spectrum  function  representation  is 
based  on  the  B66  data,  as  well  as  a  number  of  aircraft  and  towe^  measure¬ 
ments.  Tower  measurements  have  been  used  because  they  may  provide 
an  indication  of  the  behavior  of  the  longer  wavelength  spectral  components, 
as  shown  in  Figure  1 . 

The  spectrum  function  presented  in  Figure  1  is  of  the  form 

<p(k)  =  2fL/(l  +  2fkL)2  (3) 

where  k  =  f/U  is  a  wave  number  (reciprocal  wavelength),  f  is  the  frequency 
in  cycles/sec,  U  is  the  true  airspeed  in  ft/sec,  and  L  is  a  scale  para¬ 
meter  in  feet  related  to  the  terrain  roughness  and  height  above  the  ground. 

Although  Equation  (3)  does  not  provide  for  a  -5/3  power  law  be¬ 
havior  at  the  smaller  wavelength  (about  100  ft),  Figure  1  indicates  the  rela¬ 
tively  small  differences  that  exist  between  the  spectrum  estimates  for  a 
-5/3  and  a  -2  power  law.  Since  it  is  difficult  to  conclusively  prove  the 
existence  of  either  of  these  power  laws  from  present  measurements.  Equa¬ 
tion  (3)  is  considered  a  reasonable  representation  for  the  spectrum  form 
in  the  low-level  regime. 
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The  scale  parameter  L  for  the  B-66  data  «ai  estimated  for  four 
terrain  and  three  height  (200,  600,  and  1000  ft)  classifications  as  shown 
in  Figure  1  of  refs  3  and  4.  In  view  of  the  relatively  small  differences 
observed  in  the  smooth  terrain  farmlands  and  woodlands  classifications 
(both  with  respect  to  L  and  subsequent  wind  speed  characteristics),  a 
single  smooth  terrain  scale  parameter  is  believed  adequate.  The  linear 
scale  parameter  with  height  relationship  is 

L  =  hQ  +  h  1^  .  (4) 

The  values  sf  and  hQ  for  each  terrain  are  shown  in  Table  1.  In  Equa¬ 
tion  (4),  L  is  assumed  constant  above  1000  ft. 

Table  1.  Values  of  L^  and  hQ  for  each  terrain. 


Terrain  Class 

h 

o 

Smooth 

2/3 

135 

Low  Mountains 

1/2 

300 

High  Mountains 

1/8 

675 

Section  4 


VERTICAL  VELOCITY  RMS  BEHAVIOR  WITH  MEAN  WIND  SPEED 


I.  Smooth  terrain  estimate s 

Smooth  terrain  vertical  velocity  rmi  estimates  as  a  function  of 
mean  wind  speed,  based  on  the  B-66  data,  are  presented  in  refs  3  and  4 
for  three  atmospheric  stability  classifications.  Although  the  B-66  wind 
and  stability  estimates  were  not  considered  highly  reliable,  the  rms  func¬ 
tions  obtained  from  the  data  were  found  to  be  in  reasonable  agreement 
with  meteorological  tower  data. 

hi  the  analysis  of  the  climatic-wind-speed  data  conducted  in  the  pre¬ 
sent  study,  only  a  relatively  few  (about  15  percent)  neutral  stability  cases 
(see  Section  6)  were  encountered.  This  suggests  a  possible  simpliication 
to  the  model  procedure  by  including  the  neutral  cases  in  the  unstable  clas¬ 
sification.  Under  this  simplification,  relations  for  the  rms  gust  velocity 
as  a  function  of  mean  wind  speed  would  be  represented  by 


and 


Ow  =  0. 1  u  (5) 

Ow=  1.65  +.08  u  (6) 


for  stable  and  unstable  conditions,  respectively;  where  u  is  the  mean  wind 

speed  in  ft/sec.  hi  accordance  with  ref  3,  the  a  functions  are  considered 

w 

to  be  independent  of  height. 

Although  the  stable  expression  (5)  is  believed  to  be  reasonable  for 
low-to-m ode  rate  wind  speed  conditions,  it  does  not  appear  reasonable  to 
have  the  stable  a  values  exceed  the  unstable  a  values  for  large  wind  speed 
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condition.  For  this  reason, and  because  only  the  high  wind  speed  condi¬ 
tions  are  generally  of  interest  in  aircraft  gust  load  problems,  it  is  pro¬ 
posed  to  use  the  unstable  expression  to  represent  the  rms- mean- wind -speed 
behavior  in  the  present  gust  load  model. 

hi  the  application  of  Equation  (6)  the  regression  constants  should  be 
increased  to  allow  for  statistical  variability  in  aw  as  discussed  in  refs  3 

awl  4.  Based  on  954  confidence  limits  for  the  o  data  used  to  obtain  the 

w 

regression  equation.  Equation  (6)  becomes 

aw  =  2.45  +  .  12  u  (7) 

or,  in  terms  of  m/s, 

Ow  =  .  75  +  .  12  u  (8) 


2.  Rough  terrain  estimates 

As  described  in  ref  3,  the  effect  of  atmospheric  stability  could  not 
be  separated  from  the  scatter  of  the  rough  terrain  aw  estimates.  The 
least  square 8  linear  regression  lines  for  the  low  and  high  mountains  ctw 
estimates  reported  in  refs  3  and  4  were  noted  to  be  approximately  equal 
for  both  rough  terrain  categories.  The  principal  difference  between  the 
two  sets  of  data  being  the  much  greater  scatter  observed  for  the  high  moun¬ 
tain  a  estimates, 
w 

In  view  of  the  closeness  of  the  two  regression  functions,  it  is  sug¬ 
gested  that  a  single  expression  be  used  for  both  rough  terrain  categories; 
namely. 


9 


(9) 


O  =  3  +  .05  u 
w 

where  a  is  in  ft /tec.  However,  because  of  the  difference  in  the  statistic 
cal  variability  of  the  estimates,  the  equations  used  in  the  gust  load  model 
are  different.  Thus,  for  95)1  confidence  limits,  the  function  for  the  low 
mountain  case  is 

a  =  1.3+  .07  u;  (10) 

w 

whereas,  the  equivalent  expression  for  the  high  mountain  terrain  classifica¬ 
tion  is 

o  =1.8+  0.1  u  (11) 

w 

where  u  is  expressed  in  m/s. 
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Section  5 

CLIMATOLOGICAL  WIND  SPEED  DISTRIBUTIONS 

1.  Previous  data 

A  description  of  the  mean  wind  speed  distribution  characteristics 
for  smooth  and  rough  terrain  conditions  is  presented  in  refs  Z-4.  The 
smooth  te  rrain  distribution  characteristics  obtained  in  ref  2  are  based  on 
radiosonde  and  pibal  observations  obtained  at  Dayton,  Ohio  from  June  1956 
to  May  1958.  For  these  data  the  following  factors  were  considered:  atmos¬ 
pheric  stability,  seasons,  time  of  day,  and  height. 

The  results  of  additional  climatological  data  analyses  for  both 
smooth  and  rough  terrain  locations  were  reported  in  refs  3  and  4.  For 
these  results  the  smooth  terrain  data  used  were  for  June  1953  through  May 
1957,  for  1000  and  1600  EST  at  Dayton  and  for  1500  CST  at  Fort  Worth, 
Texas.  The  rough  terrain  data  were  for  Medford,  Oregon,  Las  Vegas, 
Winnemucca,  and  Ely,  Nevada  and  Boise,  Idaho  for  1200  and  2400  GMT. 

The  height  range  for  the  rough  terrain  data  extended  from  the  surface  to 
5000  m  above  the  surface. 

2.  Present  data 

To  provide  a  more  extensive  climatic  and  geographical  coverage 
for  the  climatological  wind  data,  the  data  summarized  in  Table  2  were  ob¬ 
tained  from  the  National  Weather  Records  Center  (ESSA),  Ashe  lie.  North 
Carolina: 


* 

i 


11 


Table  2.  Climatological  data  summary 


STATION  PERIOD  OF  RECORD 

TIME  HEIGHT  RANGE 

(LST)  (above  surface,  meters) 

San  Diego,  Cal. 

1958-1959 

16  00  (PST) 

0-3000 

El  Paso,  Tex. 

1958-1959 

1700(MST) 

0-5000 

Denver,  Colo. 

1958-1959 

1700(MST) 

0-6000 

Oakland,  Cal. 

1958-1959 

1600(MST) 

0-3000 

Bismark,  N.  Dak. 

1958-1959 

1700(MST) 

0-4000 

Lander,  Wyo. 

1958-1959 

1700(MST) 

0-6000 

Seattle,  Wash. 

1958-1959 

16  00  (PST) 

0-3000 

Gt.  Falls,  Mont. 

1959-1959 

1700(MST) 

0-5000 

Tucson,  Ariz. 

1960-1961 

1700(MST) 

0-4000 

Sault  St.  Marie,  Mich. 

1955-1956 

0900  (CST) 

0-762 

Topeka,  Kans. 

1955-1956 

0900(CST) 

1500(CST) 

0-762 

Montgomery,  Ala. 

1955-1956 

0900(CST) 

1500(CST) 

0-762 

Norfolk,  Va. 

1955-1956 

1000(EST) 

1600(EST) 

0-762 

3.  WinJ  distributions 

As  described  in  ref  3,  a  general  wind  probability  density  distribu¬ 
tion  was  obtained  for  the  four  conditions  listed  above  by  using  dimension¬ 
less  wind  speed  (u/u)  as  the  independent  variable,  where  u  is  the  average 
mean  wind  speed  associated  with  a  particular  condition,  e.  g. ,  height  above 
the  surface,  referred  to  as  "climatic-wind-speed". 

In  ref  3,  a  small  difference  in  the  form  of  the  dimensionless  wind 
speed  density  distribution  was  noted  for  the  smooth  and  rough  terrains. 

The  present  data  indicates  a  third  variation  in  the  probability  density  form; 
Figure  2  summarizes  the  three  distribution  forms.  For  the  present  data. 
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all  rough  terrain  stations  exhibit  greater  probability  density  values  at 
large  u/u  values  than  the  previous  rough  terrain  stations.  In  addition, 
the  Norfolk  and  Montgomery  data  also  conform  to  this  third  distributional 
form. 

The  nearest  analytic  representation  found  for  the  density  distribu¬ 
tions  is  of  the  form  described  in  ref  2: 

C2 

G,-l  -G.x 
p(x)  =  GjG2  x  e 

where  x  =  u/u;  Gj  and  G^  are  arbitrary  parameters.  This  distributional 
form  was  first  used  by  Weibull  (ref  6)  to  describe  the  observed  behavior 
of  many  physical  random  variables.-  Suggested  values  of  G^  and  G^  to 
approximate  the  wind  distribution  behavior  are  presented  in  Figure  2.  It 

i 

appears  difficult,  however,  to  obtain  a  close  approximation  to  the  observed 
data  with  this  representation. 
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Section  6 


CLIMATIC  MEAN  WIND  SPEED  CHARACTERISTICS 

1.  Wind  data 

To  use  the  climatological  wind  speed  distributions  described  in 
Section  5,  a  systematic  means  of  specifying  climatic  mean  wind  speed 
values  (u)  is  necessary  for  the  continental  United  States.  To  examine 
surface  and  upper  level  wind  speed  data  for  this  purpose,  the  NWRC  data 
presented  below  were  obtained  for  the  surface  and  for  heights  ranging  up 
to  6000  m  above  the  mean  sea  level  (MSL).  The  data  were  obtained  as 
monthly  averages  of  hourly  data  at  the  surface,  and  for  1800  and  2400  GMT 
for  the  upper  levels. 

A  summary  of  the  climatic- .vind-speed  data,  averaged  for  the  years 
listed  in  Table  2,  is  presented  in  Table  3.  Small  year-to-year  variations 
(less  than  10  percent)  were  observed  in  the  wind  speed  averages. 

2.  Rough  terrain 

(a)  Height  variation  -  In  considering  the  height  profiles  of  the  clima¬ 
tic  wind  speeds,  the  u  values  of  Table  3  were  averaged  for  the  two  hourly 
periods  shown,  since  the  differences  in  the  u  values  for  the  two  periods  are 
quite  small  at  the  upper  levels. 

In  ref  3,  the  rough  terrain  u  values  were  divided  by  the  u  surface 
values  in  an  effort  to  generalize  the  wind  speed  behavior  with  height.  Based 
on  the  additional  data  analyzed  in  the  present  study,  this  procedure  appears 
less  desirable  because  of  the  variability  observed  in  the  surface  u  values. 
Figures  3  and  4  show  the  u  values  for  14  rough  terrain  stations  and 
4  smooth  terrain  stations.  It  is  observed  that  the  variations  from  the 
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Table  3.  Climatic-wind- speed  summary, 


Height(rn) 


Station 

(LST) 

sfc 

1000 

2000 

3000 

4000 

5000 

6000 

New  York 

1900 

5.7 

10.2 

12.1 

14.6 

17.3 

20.2 

23.2 

(N.Y.) 

1300 

6.1 

9.5 

12.4 

14.2 

17.9 

20.6 

23.7 

Dayton 

1900 

4.5 

9.7 

11.0 

13.3 

16.0 

18.4 

21.0 

(Ohio) 

1300 

5.3 

10.1 

10.5 

12.8 

15.2 

17.7 

20.3 

Montgomery 

1800 

3.3 

7.0 

8.8 

10.0 

13.2 

15.2 

17.6 

(Alabama) 

1200 

4.2 

7.0 

9.4 

11.0 

13.0 

14.8 

17.4 

Topeka 

1800 

6.3 

9.2 

10.7 

12.6 

14.7 

16.8 

18.8 

(Kansas) 

1200 

6.5 

9.0 

11.2 

12.8 

14.5 

16.6 

18.6 

Great  'Falls 

1700 

6.4 

8.9 

10.8 

13.6 

16.4 

18.8 

— 

(Montana) 

1100 

6.  i 

8.8 

10.7 

13.2 

15.8 

16.5 

— - 

Long  Beach 

1600 

4.8 

— 

— 

— 

— 

— 

— 

(California) 

1000 

2.4 

— ■ 

— — 

— 

— 

— 

— 

Santa  Monica 

1600 

— 

4.1 

5.9 

8.7 

10.8 

12.7 

14.7 

(California) 

1000 

— 

4.0 

6.3 

8.8 

11.0 

12.6 

14.6 

Oakland 

1600 

5.7 

5.8 

7.2 

9.0 

— — 

— — 

— 

Station 

(LST) 

sfc 

900 

1400 

2400 

3400 

4600 

Denver 

1700 

4.8 

6.2 

6.8 

9.5 

13.0 

15.8 

(Colorado) 

1100 

3.8 

5.7 

6.9 

9.9 

13.0 

15.4 

Station 

(LST) 

sfc 

900 

1300 

1800 

2200 

3800 

El  Paso 

1700 

5.8 

6.  C 

6.8 

7.9 

10.2 

11.8 

(N.  Mexico) 

Station 

(LST) 

sfc 

1Q00 

1500 

2000 

2500 

3500 

Bismarck 

1700 

6.5 

8.8 

9.1 

10.3 

11.7 

13.8 

(N.  Dakota) 

Station 

(LST) 

sfc 

800 

1300 

2300 

3300 

4300 

Lander 

1700 

3.8 

5.2 

6. 4 

10.1 

13.2 

15.8 

(Wyoming) 

Station 

(LST) 

sfc 

700 

1200 

1700 

2200 

3200 

Tucson 

1700 

5.2 

5.7 

5.9 

6.4 

7.1 

9.7 

(Arizona) 

All  wind  speeds  in  m/s. 
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avenge  profiles  shown  are  less  than  ±  ZO  percent,  and  that  the  standard 
error  of  estimate  for  the  rough  terrain  data  is  about  one  (l)m/s. 

Of  the  14  rough  terrain  stations  analysed,  only  the  West  Coast 
Santa  Monica  and  San  Diego  profiles  characteristically  depart  from  the 
average,  suggesting  that  the  u  profiles,  in  general,  do  not  depend  on  geo¬ 
graphical  location.  It  is  interesting  to  note  that  the  rough  terrain  and 
smooth  terrain  profiles  converge  at  the  higher  elevations,  as  shown  in 
Figure  5.  ft  would  appear  that  the  obstructing  effect  of  mountainous  ter¬ 
rain  produces  a  steeper  wind  profile  for  the  first  1500  meters,  or  so, 
above  the  surface,  and  that  orographic  effect  gradually  diminishes  with 
height  and  then  disappears  at  about  the  3000  to  4000  meter  level. 

(b)  Seasonal  variations  -  To  determine  the  seasonal  variation,  the 
monthly  averages  for  six  of  the  stations*  of  Table  Z  were  divided  by  the 
corresponding  "period  of  record"  averages  listed  in  Table  3,  and  this 
dimensionless  value  was  platted  for  each  month  and  for  each  elevation 
above  the  surface  from  1000  to  5000  meters.  Since  differences  between 
the  forenoon  and  afternoon  values  were  minor,  the  two  daily  readings  were 
combined  for  each  of  the  six  stations  and  the  monthly-to-yearly  wind  speed 

ratios  (u_/u  )  plotted  for  each  height.  The  results  for  the  Z000m  eleva- 

m  an 

tion  are  shown  in  Figure  6.  The  scatter  about  the  mean  of  the  u  /u 

m  an 

values  for  the  ZOOOm  level  is  similar  to  the  scatter  observed  for  the  four 
other  elevations.  The  mean  values  for  all  five  elevations  are  shown  in 


The  Denver  station  elevation  of  1600m(MSL)  precluded  conveniently  in¬ 
cluding  this  station  in  the  even  numbered  height-above-station  classifica¬ 
tion. 
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Figure  8,  where  the  four  single  year  distributions  are  shown  for  the  New 
York  station  (1961  to  1964)  for  the  2000  m  level. 

3.  Smooth  terrain 

(a)  variarinn  .  The  monthly  variations  of  the  um/uan  wiDd 
speed  ratios  at  the  surface  are  shown  in  Figure  9  for  the  four  smooth  ter¬ 
rain  stations.  The  year-to-year  variations  for  the  New  York  station  are 
shown  in  Figure  10.  Although  the  year-to-year  single  station  variations 

at  the  surface  are  seen  to  be  greater  than  the  upper  level  year-to-year  dis¬ 
tribution  (Figure  8),  the  station-to-staticm  variability  of  the  averaged 
utTi/uaw  values  of  Figure  9  are  comparable  to  those  of  Figure  6  for  the 
upper  levels.  The  principal  difference  between  the  surface  and  upper  level 
distributions  is  the  greater  skewness  of  the  surface  distribution.  The  sea¬ 
sonal  variation  of  the  wind  speeds  at  the  300  m  level  are  described  under 
(c)  below. 

(b)  Diurnal  variation  -  Surface  wind  speeds  obtained  from  hourly 
data  are  shown  in  Figures  11  and  12  for  four  smooth  terrain  stations  and 
Long  Beach.  These  figures  indicate  that  the  year-to-year  variations  for 
all  stations  are  comparatively  small.  In  Figure  13  the  envelope  of  the 
monthly  averages  are  shown  for  Dayton  for  1961. 

To  generalize  the  diurnal  wind  speed  behavior,  the  average  hourly 
distributions  of  Figures  11  and  12  were  divided  by  the  1500  local  standard 
time  wind  speeds,  which  represent  the  approximate  daily  maximum  values. 
The  results  of  this  normalization,  shown  in  Figure  14,  reveals  a  fairly  con¬ 
sistent  pattern  of  diurnal  variations  from  about  0800  to  1700  (LST),  although 
both  the  early  morning  and  evening  wind  speed  averages  are  seen  to  differ 

by  as  much  as  10$  from  the  mean  at  particular  stations. 
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The  drastic  change  in  the  behavior  of  the  Long  Beach  surface  winds 
is  believed  to  be  produced  by  the  proximity  of  the  ocean  and  the  California 
current  off  the  coast  of  southern  California.  The  over-water  trajectories 
of  the  westerly  winds  provide  relatively  small  wind  speeds  through  most 
of  the  day;  however,  the  development  of  an  afternoon  sea  breese  in  this  re¬ 
gion  may  be  the  reason  for  the  increase  in  the  mean  wind  speed  values  dur¬ 
ing  the  mid -afternoon  period,  raising  the  wind  speed  values  to  about  the 
average  of  the  other  continental  stations,  hi  Table  4  the  1500(LST)  average 
wind  speeds  are  listed  for  both  rough  and  smooth  terrain  stations;  the  aver¬ 
age  for  the  smooth  terrain  stations  is  just  under  5. 5  m/s,  while  the  rough 
terrain  station  average  is  6  m/s.  The  Long  Beach-San  Diego  average  is 
5  m/s.  As  a  matter  of  interest,  it  was  ascertained  that  the  diurnal  surface 
wind  variations  for  the  rough  terrain  stations  were  substantially  similar  to 
the  smooth  terrain  variations. 

(c)  Atmospheric  stability  conditions  -  hi  ref  2,  the  effect  of  atmos¬ 
pheric  stability  on  the  mean  wind  speed  values  was  examined  in  consider¬ 
able  detail  for  two  smooth  terrain  locations:  Dayton,  Ohio  and  Fort  Worth, 
Texas,  hi  the  present  study,  stability  wind  speed  data  for  three  additional 
smooth  terrain  stations  were  obtained.  The  stability  criteria  used  for 
these  data  was  the  temperature  difference  between  the  300m  level  and  sur¬ 
face  station.  A  temperature  difference  less  than  5*F  was  considered  stable 
between  5*F  and  6*F  neutral,  and  greater  than  6*F  unstable.  Soundings 
were  available  for  both  morning  and  afternoon  at  the  three  stations.  A 
summary  of  the  frequency  of  occurrence,  mean  wind  speed,  and  stability 
distribution  (in  percent)  for  each  sounding  period  is  summarised  in  Table  5 


Table  4.  Summary  of  climatic -wind-speed  values  at  the  surface  for 
1500  LST. 


Station 

Wind  Speed 
(m/s) 

Station 

Wind  Speed 
(m/s) 

Dayton 

(Ohio) 

5.5 

El  Paso 
(Texas) 

6.2 

Norfolk 

(Virginia) 

5.2 

Tucson 

(Arizona) 

5.7 

Topeka 

(Kansas) 

6.3 

Lander 

(Wyoming) 

4.9 

Montgomery 

(Alabama) 

4.4 

Great  Falls 
(Montana) 

7.0 

New  York 
(New  York) 

6.3 

Denver 

(Colorado) 

5.0 

Long  Beach 
(California) 

4.8 

Bismarck 
(N.  Dakota) 

7.0 

San  Diego 

5.2 

Oakland 

6.2 

(California) 

(California) 

Sault  St.  Marie 
(Michigan) 

5.1 

Seattle 

(Washington) 

5.8 

Table  5.  Smooth  terrain  climatic-wind-speed  ratios  for  three  stability 


conditions 

a 

u 

300^usfc 

Frequency 

“sfc 

(percent) 

Station 

(LST) 

_S 

N 

_U 

S 

N 

_U 

_S 

N 

U 

Dayton 

0930 

2.1 

1.7 

1.7 

69 

11 

20 

3.5 

4.4 

3.7 

(Ohio) 

1530 

1.8 

1.5 

1.5 

38 

26 

36 

4.5 

3.8 

4.2 

Norfolk 

1000 

1.6 

1.5 

1.5 

65 

10 

25 

5.1 

5.5 

5.8 

(Va.) 

1600 

1.6 

1.4 

1.3 

57 

16 

26 

5.3 

5.4 

5.0 

Montg. 

0900 

2.1 

1.7 

1.5 

50 

12 

38 

3.2 

3.7 

3.2 

(Ala.) 

1500 

1.8 

1.4 

1.5 

29 

17 

54 

4.2 

4.3 

3.6 

Topeka 

0900 

1.6 

1.5 

1.4 

73 

8 

19 

5.1 

5.0 

5.5 

(Kane. ) 

1500 

1.3 

1.2 

1.2 

30 

22 

48 

6.1 

6.7 

6.3 
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for  the  five  smooth  terrain  stations,  la  addition.  Table  5  presents  300 in¬ 
to- surface  wind  speed  ratios  for  each  stability  condition. 

In  Table  5,  the  wind  speed  ratios ( u3oo^usfc^ *or  ****  neutra^  an^  un~ 
stable  conditions  are  seen  to  vary  much  less  than  the  same  wind  speed  ratios 
for  the  stable  classification.  Moreover,  if  the  wind  speed  ratio  is  plotted 
with  the  surface  wind  speed,  as  in  Figure  17,  the  stable  ratios  are  shown 
to  decrease  with  increasing  ug^c  values.  The  wind  speed  ratios  for  both 
the  neutral  and  unstable  conditions,  on  the  other  hand,  do  not  correlate  as 
well  with  the  values,  h  fact,  below  about  6m/s  die  neutral  and  un¬ 

stable  ratios  vary  only  about  10  percent  from  a  mean  value  of  about  1.5. 

The  relatively  small  number  of  neutrai  stability  cases  indicated  by 
the  data  of  Table  5,  and  the  small  difference  in  the  u300^u8fc  ratios  for 
the  neutral  and  unstable  cases,  suggest  the  possibility  of  combining  these 
two  stability  classes  into  one  grouping  for  modeling  purposes.  On  this 
basis,  the  diurnal  distribution  of  the  relative  number  of  unstable  (neutral 
and  unstable  combined)  cases  are  shown  in  Figure  18  as  a  function  of  local 
standard  time.  Although  these  results  do  indicate  variations  from  station 
to  station,  no  specific  geographical  trend  is  discernible.  The  approximate 
behavior  indicated  in  the  figure  is  suggested  to  represent  the  relative  oc¬ 
currence  of  unstable  stability  conditions  between  the  hours  of  0600  and 
1800  LST. 

The  wind  speed  ratio  values  in  Figure  17  (stable  cases  only)  can  be 

=  2.9  -  .  26  usfc  (12) 


approximated  by  the  expression 


sfc 
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For  Intermediate  altitudes  of  150  and  60  meters,  the  profile  characteris¬ 
tics  reported  in  refs  2  and  3  are  used  to  approximate  ui5()/u8fc  anc* 
u60^usfc  ra^os*  Using  Figure  15  of  ref  2,  and  the  data  of  Figure  17,  the 
wind  speed  ratios  for  the  stable  condition  at  150  and  60  meters  are  repre¬ 
sented  by 

=2.3  -.18  usfc  .  (13) 

usfc 

and 

=  1.8  -.10us£c  .  <14) 

usfc 

Under  neutral  to  unstable  stability  conditions  the  average  wind  speeds 
do  not  vary  appreciably  between  the  500  ft  (150m)  and  1000ft  (300m)  levels, 
as  shown  by  Figure  17  of  ref  2.  At  150  m  the  ui5o^uflfc  rat*°  from  ref  2 
is  about  1.5,  in  approximate  agreement  with  the  data  of  Table  5.  Hence, 
these  values  are  used  in  the  present  model. 

The  seasonal  variation  of  the  climatic-wind-speeds  at  300  m  for 
stable  and  unstable  conditions  were  considered  for  the  available  morning 
and  afternoon  time  periods  of  1500  and  2100  GMT.  For  the  stable  data, 
the  monthly-to-yearly  wind  ratios  did  not  vary  significantly  for  the  two  time 
periods;  however,  the  seasonal  distributions  for  the  unstable  data  are  shown 
in  Figure  16  to  be  decidedly  dissimilar  for  the  two  time  periods.  The  wind 
speed  ratios  for  the  stable  condition  are  shown  in  Figure  15.  The  smooth 
curves  in  these  figures  are  suggested  representions  for  the  three  set  of  data. 
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Section  7 


GENERALIZED  LOAD  EXCEEDANCE  CURVES 


Based  on  the  results  of  Sections  4-6,  generalized  load  exceedance 
design  curves  can  be  obtained  from  Equation  (1).  For  this  purpose.  Equa¬ 
tion  (1)  is  expressed  as 

5  £(ffw>  exp  ffl}  daw  (15) 

o  2<t 

e  w 

where  x  =  Y/A. 

For  specified  values  of  Equation  (15)  is  evaluated  for  an  aver¬ 
age  probability  density  function  of  mean  wind  speed  (Figure  2)  and  a  range 
of  values  of  a,  b,  and  u,  representing  the  cr^  dependence  on  mean  wind 
speed: 

a  =  a  +  b  u 
w 

The  a  and  b  values  were  described  in  Section  4  for  the  smooth  and  rough 
terrain  conditions;  the  characteristics  of  u  were  described  in  Section  6. 
Since  both  b  and  u  vary,  the  product  bu  is  represented  by  fJ  below. 

Solutions  for  Equation  (15)  were  obtained  for  four  specific  values  of 
X:  7,  14,  21,  and  28.  These  solutions  were  obtained  for  (3  values  ranging 
from  0.  5  to  1.8  for  parametric  values  of  a  of  0,  0.5,  1.0,  1.5,  and  2.0. 
These  solutions  are  shown  in  Figure  20. 

To  obtain  solutions  for  values  of  x  and  a  intermediate  to  those 
shown  in  Figure  20,  cross-plots  can  be  constructed.  For  instance,  to 
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evaluate  intermediate  values  of  a  (for  a  given  N/NQ  values  are  plotted 
as  a  function  of  a  for  a  given  value  of  fi;  intermediate  values  of  a  are  then 
obtained  from  the  curve  drawn  through  the  five  a  values  provided  by  Figure 
20.  Similarly,  intermediate  values  of  Y/A  can  be  obtained  from  a  plot  of 
N/Nq  versus  \  for  given  values  of  a  and  0. 
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Section  8 


MODEL  PROCEDURES 


In  Sections  3  and  4,  the  spectrum  shape  function  and  vertical  velo¬ 
city  rms  characteristics  with  respect  to  terrain  roughness  were  described. 
In  Sections  5  and  6  (and  refs  2  and  3),  the  climatological  wind  speed  prob¬ 
ability  distributions  and  climatic-wind-speed  characteristics  with  respect 
to  such  variables  as  height,  seasons,  atmospheric  stability,  time  of  day, 
and  terrain  roughness  were  presented.  In  Section  7,  generalized  load  ex¬ 
ceedance  curves  were  presented  to  permit  relatively  easy  evaluation  of 
N/Nq  as  a  function  of  Y/A  for  the  parameters  (a,  0)  representing  the 
dependence  on  terrain,  height  and  climatic-wind-speed  conditions.  In  this 
section,  an  outline  of  procedures  is  presented  for  utilizing  detailed  clima¬ 
tological  wind  speed  characteristics  for  estimating  gust  load  exceedance 
probabilities.  Also,  an  illustration  of  the  application  of  these  procedures 
is  presented  for  the  rough  terrain  case. 

1.  Smooth  terrain  procedure 

The  average  number  of  peak  values  per  unit  distance  exceeding  a 

given  load  level  is  given  by  Equation  (15).  For  a  given  height  above  the 

2 

surface  and  a  given  true  airspeed,  A  and  Nq  depend  on  the  normalized 
spectrum  function  and  the  aircraft  dynamic  response  function  T(ifl).  For 
a  particular  height  and  aircraft  flight  condition  (true  airspeed,  weight,  etc), 
the  procedure  for  evaluating  Equation  (15)  for  smooth  terrain  conditions  is 
described  below,  and  in  flow  chart  form  in  Figure  21.  For  the  present  it 
is  assumed  that  all  smooth  terrain  conditions  are  sufficiently  similar  to 
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permit  using  the  statistical  averages  obtained  from  the  present  study.  It 
is  possible,  however,  that  for  some  applications  it  maybe  desirable  to 
evaluate  the  model  for  specified  smooth  terrain  locations. 

1)  The  "operational  day"  is  divided  into  a  number  of  "diurnal  time  units", 
disregarding,  for  the  present,  the  relative  amount  of  time  spend  in  each 
unit.  Selection  of  time  units  or  blocks  is  based  on  Figure  14. 

2)  The  daytime  maximum  surface  climatic-wind-speed  value  for  a  specific 
smooth  terrain  location  is  selected  (see  Table  4),  or  a  value  of  5.  5  m/s 
is  used  to  represent  general  smooth  ter..  ;in  conditions.  Wind  speed 
values  corresponding  to  diurnal  time  units  of  step  (1)  are  obtained  from 
Figure  14. 

3)  Select  "seasonal  blocks"  on  the  basis  of  Figures  15  and  16  to  represent 
seasonal  changes  in  climatic- wind -speed.  For  each  stability  condition 
multiply  the  wind  estimates  of  step  (2)  by  the  appropriate  values  from 
Figures  15  and  16.  (The  monthly  changes  in  the  wind  speeds  at  300m, 
represented  by  these  figures  for  stable  and  unstable  conditions,  are  also 
used  to  represent  the  changes  at  the  15C  and  60  meter  levels). 

4)  For  stable  atmospheric  stability  conditions,  the  climatic-wind-speed  values 
at  a  given  flight  level  are  obtained  from  the  surface  values  of  step  (3)  and 
one  of  the  Equations  (12-14)  based  on  Figure  17.  For  unstable  stability 
conditions,  the  flight  level  climatic-wind-speeds  are  obtained  by  increas¬ 
ing  the  surface  wind  speeds  by  50  percent  for  the  300m  and  150  m  levels, 
and  by  35  percent  for  the  60  m  level,  in  accordance  with  Section  6. 
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5)  Using  the  a  values  of  step  (4)  and  the  vertical  velocity  rms  function 
[Equation  (7)  and  (8)],  estimate  a  and  0  for  the  stable  and  unstable 
conditions. 

6)  For  selected  Y/A  values,  £<  Ermine  N/Nq  for  a  and  0{step  (5)]  from 
Figure  20. 

7)  Figure  18  is  used  to  divide  each  diurnal  time  unit  of  step  (1)  into  a 
stable  portion;  and  a  relative  exposure  weight  (REW)  is  assigned  to  re> 
present  the  percent  of  smooth  terrain  flight  time  at  each  diurnal  time 
unit  for  each  location,  height,  and  season. 

8)  Using  the  numerical  weights  determined  in  step  (7),  N/Nq  values  are 
determined  for  stable  and  unstable  conditions.  (N^  ;N^  in  Figure  21). 

9)  Add  the  N/N  values  of  step  (8) :  (N.  +  N,)/N  =  N./N  . 

10)  Repeat  steps  (1)  through  (9)  for  each  time  and  seasonal  unit  selected 
(for  a  given  height  and  location): 

L  Nij<Y)  =  N-  *WY) 

o  ..  o 

where  i  and  j  represent  time  and  seasonal  units  and  z  and  k  height 
and  location. 

11)  Repeat  steps  (1)  through  (10)  for  each  location  desired.  Thus, 

N/No<Y>  =  N  l  Nzk(Y>  =  N  Nz<Y» 
o  ,  o 

k 

where  Nz(Y)  represents  the  summation  over  j  and  k  for  a  single  height 
z. 

12)  Evaluate  Nq  for  height  above  surface  used  in  steps  (1)  through  (11). 
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13)  Repeat  steps  (1)  through  (11)  for  each  additional  height  desired.  Sam 
the  N(Y)  resalts  for  all  heights  selected: 

^^Smooth  terrain  *  I  N.'T'  • 

a 

2.  Rough  terrain  procedure 

The  procedure  followed  for  the  rough  terrain  cases  is  similar  to 
the  smooth  terrain,  except  that  neither  atmospheric  stability  nor  diurnal 
variations  of  the  upper  level  wind  speeds  need  be  considered.  For  a  parti¬ 
cular  height  and  aircraft  flight  condition  the  rough  terrain  procedure  is  des¬ 
cribed  below  and  in  flow  chart  form  in  Figure  22.  As  in  the  smooth  terrain 
case,  the  climatic-wind-speed  statistics  are  not  related  to  part*~ular  loca¬ 
tions;  although,  as  for  smooth  terrains,  it  is  possible  to  consider  particular 
terrain  locations,  when  necessary. 

1)  Apportion  flight  time  over  rough  terrain  between  low  and  high  mountain 
conditions;  and  determine  the  average  flight  levels  (above  the  surface 
stations)  for  the  terrain  locations  selected.  (See  discussion  below). 

2)  Use  Figure  3  to  estimate  climatic-wind-speeds  for  the  heights  selected 
in  step  (1). 

3)  Use  Figure  7  to  adjust  the  wind  speed  values  obtained  from  step  (2)  for 
seasonal  variations. 

4)  Using  the  u  values  of  step  (3)  and  the  rough  terrain  rms  functions  [Equa¬ 
tion  (10)  and  (11)],  estimate  a  and  p  for  each  function. 

5)  Determine  values  of  A  for  heights  z.  selected  for  high  and  low  mountain. 

6)  For  values  of  A  of  step  (5)  and  a  selected  Y,  determine  N/N  for  a 

o 

P  of  step  (4)  from  Figure  20. 
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7)  Assign  relative  exposure  weights  for  each  height  to  represent  percent 
of  rough  terrain  flight  time  at  height  z^  above  the  surface.  Designate 
this  Cj-  Also  assign  a  relative  weight  for  the  percent  of  rough  terrain 
flight  time  over  high  mountain  terrain.  Designate  this  r;  (1  -  r)  is  then 
the  percent  of  rough  terrain  flight  time  over  low  mountain  terrains. 

8)  Using  the  relative  weights  of  step  (7),  the  low  mountain  and  high  moun¬ 
tain  N/Nq  values  are  obtained  (N^/Nq  and  N^/Nq  in  Figure  22). 

9)  Evaluate  Nq  for  the  heights  selected  for  the  low  and  high  mountain  (z.). 

10)  Multiply  results  of  step  (8)  by  the  Nq  values  of  step  (9). 

11)  Add  the  load  exceedance  probabilities  Nj  and  for  the  low  and  high 
mountain  cases. 

12)  Repeat  steps  (1)  through  (11)  for  each  additional  height  and  seasonal 
block  desired. 

Nk  ■  I 

13)  Repeat  steps  (1)  through  (12)  for  specific  locations  desired. 

N(Y»Rough  terrain  =  i  Nk<Y> 

k 


Finally,  to  combine  the  load  exceedance  estimates  N(Y)gmoo^  and 
N(Y)Rougk  ,  it  is  necessary  to  provide  relative  weights  for  flight  time 
divided  between  smooth  and  rough  terrains.  Thus, 


NiYl^nNIY^ooth*  <l-n)N(Y)Rough 

where  n  is  the  percent  flight  time  over  smooth  terrains. 
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3.  Illustration  of  application  of  model 

As  an  illustration  of  detailed  application  of  the  gust  load  model  sum¬ 
marised  in  Figures  21  and  22,  the  following  rough  terrain  conditions  will 
be  assumed: 

Zj  (low  mount)  =  1000  m 

z_  (high  mount)  =  3000  m 

m  (seasonal  unit)  =  January 

Y/A  (exceedance  value  desired)  =  14 

gj  (percent  flight  time  at  z):  gj  =  .4;  g2  =  .2 

r  (percent  flight  time  over  high  mountains)  =  .4 

Nq  (characteristic  frequency)  =  .  5 

Using  the  above  data  and  the  flow  chart  of  Figure  22,  the  solution 
for  N^(Y)  and  N^(Y)  are  easily  obtained.  The  steps  performed  for  this 
solution  are  indicated  in  Figure  23,  where  the  steps  indicated  in  Figure  22 
have  been  evaluated  in  accordance  with  the  assumed  data  above. 

4.  Rough  terrain  classification 

To  determine  the  height  above  the  surface  to  use  in  the  model  it  is 
necessary  to  consider  flight  path  profiles  in  relation  to  rough  terrain  pro¬ 
files.  For  this  task  some  criteria  to  describe  the  difference  between  the 
low  and  high  mountain  terrain  classification  is  necessary.  Inasmuch  as 
this  terrain  differentiation  was  originally  used  in  the  B-66  program,  the 
rough  terrain  track  profiles  reported  by  Douglas  are  shown  in  Figure  12. 

bt  examining  these  profiles,  in  conjunction  with  rms  vertical  velo¬ 
city  data,  it  appears  that  the  low  mountain  (lass  roughness)  classification 
might  include  profiles  21  and  22  at  Kirtland  as  limiting  profiles  for  the  low 
mountain  category;  tracks  11  and  12  at  Kirtland  and  11  at  Edwards  appear 
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to  be  very  nearly  equal  in  roughness  to  21  and  22  for  Kirtland.  The  re¬ 
maining  high  mountain  prefiles  (21  Edwards  and  Shaw,  23  Edwards  and 
Kirtland,  and  24  Edwards)  provide  considerably  greater  scatter  in  the  rms 
gust  velocity  data.  Based  on  these  rough  terrain  profiles  it  is  suggested 
that  the  low  mountain  terrain  classification  include  terrain  profiles  whose 
peak-to- valle y  differences  (over  20  to  30  mile  segments)  lie  between  500 
to  2000  ft;  and  that  the  high  mountain  classification  include  all  profiles  with 
peak-to- valley  differences  exceeding  2000  ft. 

For  estimating  rough  terrain  topographical  profiles,  two  series  of 
TTSAF  Navigational  Charts  are  available.  These  are  (1)  USAF  Jet  Naviga¬ 
tion  Charts,  and  (2)  USAF  Operational  Navigation  Charts.  They  are  avail¬ 
able  from  the  Department  of  the  Air  Force,  Headquarters,  Aeronautical 
Chart  and  Information  Center,  Second  and  Arsenal,  St.  Louis,  Missouri 
63118. 


Section  9 


PRELIMINARY  COMPARISON  WITH  B-52  DATA 

As  a  preliminary  comparison  of  the  turbulence  model  with  avail* 
able  aircraft  data,  load  history  exceedance  values  reported  by  Boeing  for 
the  B-52  aircraft  were  considered. 

Although  it  was  not  possible  in  the  present  program  to  conduct  a 
detailed  comparison  of  the  Boeing  load  data  with  the  turbulence  model, 
rough  comparisons  were  made  based  on  N(Y)/No  versus  Y/A  summaries 
presented  in  refs  9-11. 

Before  discussing  the  results  of  this  comparison,  some  of  the  prob¬ 
lems  associated  with  the  B-52  data  for  the  G/H  and  B-F  aircraft  are  des¬ 
cribed.  For  present  purposes  there  are  two  principal  limitations  that  re¬ 
duce  the  usefulness  of  the  B-52  data  as  a  test  for  the  turbulence  model: 

1)  The  oil  burner  routes  were  flown  at  a  minimum  terrain-clearance 
altitude  of  800ft  over  the  plains,  and  several  thousand  feet  over  the  moun¬ 
tains.  The  poker  deck  routes  were  flown  down  to  500ft  terrain  clearance. 
The  poker  deck  data,  however,  represents  only  about  7  percent  of  the  total 
low-level  experience  and  detailed  summaries  of  these  flights  are  not  avail¬ 
able. 

2)  The  B-52  was  operating  under  flight  restrictions  which  precluded 
recording  some  of  the  more  severe  gust  load  data. 

Keeping  these  facts  in  mind,  summaries  of  N(Y)/Nq  versus  Y/A 
were  obtained  from  ref  9-11.  The  summaries  are  shown  in  Figure  24. 

The  following  is  observed  about  the  B-52  data:  1)  The  average  terrain-clear- 
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ance  heigfcft  for  the  low-level  data  is  about  2500  ft;  2)  For  the  available 
800  ft  data,  the  800  ft  poker  deck  exceedances  are  slightly  greater  than 
the  800  ft  oil  burner  exceedances  -  apparently  reflecting  the  influence  of 
contour  flying  in  the  poker  deck  data;  3)  The  800  ft  exceedances  are  con¬ 
siderably  greater  than  the  2500ft  data,  indicating  the  influence  of  the 
earth's  surface  on  the  turbulence  intensities;  4)  The  relatively  high  ter¬ 
rain-clearance  altitude  flown  for  the  ail  burner  routes  appears  to  have 
washed  out  the  terrain  influence,  as  shown  by  the  very  close  agreement 
between  the  Plains  and  Mountains  data  report  in  ref  11 . 

hi  order  to  obtain  a  rough  comparison  at  this  time  of  the  NYU 
model  with  the  B-52  data  no  attempt  was  made  to  take  into  account  detailed 
flight  conditions.  Instead,  it  is  assumed  that  on  the  basis  of  the  u  clima¬ 
tic  wind  speed  behavior  described  in  Section  6  an  average  u  value  between 
7  and  8m/s  can  be  taken  as  representative  of  wind  speed  conditions  en¬ 
countered  in  the  B-52  program  for  the  800  ft  data.  Note  that  the  Y/A  values 
for  the  model  are  in  m/s  and  must  be  multiplied  by  3.28  to  correspond  to 
Y/A  of  the  measured  load  data.  The  model  N(Y)/Nq  estimates  for  Y/A 
values  of  23  and  46  are  indicated  in  Figure  24  for  an  assumed  smooth  ter¬ 
rain  u  value  of  7.  5m/s  at  an  altitude  of  800  ft  above  the  surface. 
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Section  10 


DISCUSSION  OF  MODEL  APPLICATION 

An  outline  of  procedures  for  applying  the  low  altitude  turbulence 
model  to  estimate  aircraft  gust  loads  was  presented  in  Section  8.  The 
actual  experience  of  using  the  model  will  determine  the  necessity  of  modi¬ 
fying  or  simplifying  the  steps  outlined.  A  very  rough  comparison  of  the 
model  was  made  in  Section  9  with  a  limited  amount  of  B-52  load  exceedance 
data  obtained  over  smooth  terrain  at  an  approximate  terrain-clearance  alti¬ 
tude  of  800ft.  In  general,  two  kinds  of  "model  experience"  exercises 
should  be  considered. 

First,  the  sensitivity  of  the  model  should  be  tested  to  evaluate  the 
effect  of  (1)  modifying  the  climatic-wind-speed  probability  density  function 
for  smooth  terrain  and  rough  terrain  conditions, and  (2)  changing  the  number 
of  incremental  divisions  used  to  represent  the  physical  conditions  of  height 
and  diurnal/seasonal  blocks. 

Second,  the  operational  histories  of  low-altitude  gust  load  programs, 
such  as  the  B-52  flight  loads  program,  should  be  used  as  input  data  to  the 
turbulence  model  to  provide  such  details  as  airspeeds,  aircraft  weights, 
terrain-clearance  altitudes,  track  locations,  time  of  day  and  season.  The 
predicted  model  results  can  then  be  compared  directly  with  the  measured 
service  load  data.  Such  comparisons  would  serve  also  to  provide  statisti¬ 
cal  criteria  for  the  confidence  limits  required  to  describe  the  rms  gust 
velocity  as  a  function  of  mean  wind  speed. 
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Section  11 

SUMMARY  AND  RECOMMENDATIONS 

In  the  present  study  the  following  results  were  obtained  in  the 
development  of  a  low-level  turbulence  model: 

(1 )  Probability  density  distribution  functions  of  mean  wind  speed  to  clima¬ 
tic  mean  wind  speed  (u/u)  were  determined  for  smooth  and  rough  ter¬ 
rain  conditions. 

(2)  Climatic  mean  wind  speed  variations  with  height  were  described  for 
smooth  and  rough  terrain  conditions. 

(3)  Monthly  variations  of  the  climatic  mean  wind  speed  ratios  (u  /u  ) 

m  an 

were  obtained  for  smooth  and  rough  terrain  conditions.  For  smooth 
terrains,  the  monthly  variations  were  obtained  for  the  surface  and  300 
meter  levels.  For  rough  terrains,  the  upper  level  wind  speed  ratios 
were  found  essentially  independent  of  height. 

(4)  The  diurnal  variation  in  the  relative  number  (percent)  of  unstable  at¬ 
mospheric  stability  conditions  occurring  for  smooth  terrain  stations 
was  estimated.  In  this  estimate  the  neutral  conditions  were  included 
in  the  unstable  classification,  since  only  about  15  percent  of  the  cases 
(on  the  average)  were  found  in  the  neutral  category. 

(5)  For  smooth  terrain  conditions,  the  daytime  variations  in  the  climatic 

mean  wind  speeds  were  described  and  a  representative  function 

(u„/u _ )  was  obtained.  The  maximum  climatic  mean  wind  speed 

n  max 

values  were  found  to  occur  at  about  1500  LST. 

(6)  Based  on  the  climatic  mean  wind  speed  characteristics  obtained  in  the 

study,  procedures  were  outlined  for  applying  the  low-level  turbulence 

model  to  gust  load  exceedance  estimates. 
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(7)  A  rough  comparison  of  the  gust  load  exceedance  values  obtained  by 
B-52  aircraft  operating  at  approximate  800  ft  terrain-clearance  heights 
over  plains  and  the  turbulence  model  values  for  smooth  terrain  condi¬ 
tions  was  made. 

On  the  basis  of  the  low-level  model  developed,  the  following  recom¬ 
mendations  are  made: 

(1)  The  sensitivity  of  the  model  should  be  tested  to  determine  (a)  the 
number  of  incremental  steps  (or  blocks)  required  to  adequately 
define  the  environmental  conditions  (e.g.  height;  diurnal  and  sea¬ 
sonal  variations)  and  (b)  the  effect  of  using  separate  probability 
wind  speed  density  functions  to  represent  smooth  and  rough  terrain 
conditions. 

(2)  The  model  predictions  should  be  compared  with  available  low-alti¬ 
tude  load  exceedance  data.  For  such  comparisons  the  environmen¬ 
tal  conditions  of  time  of  day  and  year,  and  topography  associated 
with  measured  data  should  be  available. 

(3)  Additional  aircraft  measurements  of  rms  gust  velocity  and  associ¬ 
ated  wind  speed  should  be  obtained  to  more  adequately  establish 
the  relationship  between  the  two  quantities. 
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FIG.  2.  WIND  SPEED  PROBABILITY  DENSITY  DISTRIBUTIONS  FOR 
SMOOTH  AND  ROUGH  TERRAINS. 


Winnerr.ucco 
Las  Vegas 


u  m/s 


<0 


«  10  es»  —  o 

(sj»|9w  0001 )  suoijois  eoojjns  0M>qv  *mB|9H 


41 


8  9  10  II  12  I  2  3  4  5  6  7 

m  (month) 


FIG.  6.  MONTHLY-TO-YEARLY  CLIMATIC- WIND-SPEED  RATIOS 
AT  THE  2000  M  LEVEL  FOR  SIX  SURFACE  STATIONS 
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FIG,  8.  MONTHLY-TO-YEARDY  CLIMATIC-WIND-SPEED  RATIOS  FOR 
FOUR  YEARS  AT  2000  M  ABOVE  THE  NEW  YORK  STATION 
(READINGS  AT  1300  LST). 
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FIG.  9 
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.  MONTH JuY-TOYEARLY  CLIMATIC-WIND-SPEED  RATIOS  AT 
THE  SURFACE  FOR  FOUR  SMOOTH  TERRAIN  STATIONS. 
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FIG.  10.  MONTHLY-TO-YEARLY  CLIMATIC- WIND-SPEED  RATIOS  FOR 
FOUR  YEARS  NEAR  THE  SURFACE  FOR  NEW  YORK  (OBSER¬ 
VATIONS  OBTAINED  AT  1900  i£Tf. 


H  ( hour ,  LST ) 


FIG.  11.  HOURLY  CLIMATIC-WIND-SPEED  VALUES  NEAR  THE  SURFACE 
FOR  TOPEKA,  KANSAS  AND  MONTGOMERY,  ALABAMA  FOR 
1962-1964. 
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FIG.  13.  ENVELOPE  OF  HOUR1 Y-TO-MONTHLY  CLIMATIC- WIND-SPEED 
RATIOS  NEAR  THE  SURFACE  FOR  DAYTON.  OHIO  FOR  1961. 
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FIG.  14.  HOURLY-TOMAXIMUM  HOURLY  CLIMATIC-WIND-SPEED 
RATIOS  NEAR  THE  SURFACE  FOR  FIVE  STATIONS. 
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FIG.  17.  CLIMATIC- WIND-SPEED  RATIOS  uz/u8fc  FOR  STABLE  CONDITIONS 
OVER  SMOOTH  TERRAINS. 
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FIG.  18.  VARIATION  OF  UNSTABLE  (INCLUDING  NEUTRAL)  STABILITY 
CONDITIONS  FOR  DAYTIME  HOURS. 
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FIG.  20.  GENERALIZED  LOAD  EXCEEDANCE  CURVES 
(a)  Y/A  =  7;  Y/A  =  14 
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for  aircraft  gust  load  application  is  described.  The  basic  data  used  for 
this  model  are  the  B-66B  low-level  gust  data  and  climatological  wind  data 
provided  by  the  National  Weather  Records  Center,  Asheville,  N.  C.  In 
preceding  studies  the  turbulence  spectrum  functions,  rms  gust  velocity/ 
mean  wind  speed  functions,  and  mean  wind  speed  distribution  functions 
were  determined  for  variations  in  height,  atmospheric  stability,  and  ter¬ 
rain  roughness  conditions.  In  this  final  study  needed  to  complete  the  tur¬ 
bulence  model,  climatological  wind  statistics  are  used  to  relate  average 
mean  wind  speed  characteristics  to  terrain,  height,  thermal  stability,  time 
of  day  and  seasonal  variations.  The  data  represent  locations  throughout 
the  United  States.  Based  on  results  obtained  from  the  climatological  wind 
statistics,  procedures  are  outlined  for  applying  the  turbulence  model  to 
estimate  aircraft  gust  load  exceedances  for  a  specified  low  altitude  opera¬ 
tional  history.  A  preliminary  comparison  of  the  model  and  B-5Z  service 
load  data  is  made  for  an  300  ft  terrain-clearance  altitude. 
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